Identification of flow-background to subtract in jet-like azimuthal correlation 
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We derive an analytical form for flow-background to jet-like azimuthal correlation in a cluster 
approach. We argue that the elliptic flow parameter to use in jet-correlation background is that from 
two-particle method excluding non-flow correlation unrelated to the reaction plane, but including 
cross-terms between cluster correlation and cluster flow. We verify our result with Monte Carlo 
simulations. We discuss implications of our finding in the context of jet-like correlations from STAR 
and PHENIX. 
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I. INTRODUCTION 



Jet-like angular correlation studies with high transverse momentum (j>t) trigger particles have provided valuable in- 
formation on the properties of the medium created in relativistic heavy- ion collisions [H,0]- I n such studies, correlation 
functions are formed in azimuthal angle difference between an associated particle and a high px trigger particle, which 
preferentially selects (di-)jet. One important aspect of these studies is the subtraction of combinatorial background 
£SJ ■ which itself is non-uniform due to anisotropic particle distribution with respect to the reaction plane- both the trigger 
particle and the associated particles are correlated with the common reaction plane in an event. One critical part is 
1 ^ '. to determine flow parameters, mainlyelliptic flow (i>2), to use in constructing background. 

' There are many vi measurements [3, 4]. They contain various degrees of non-flow contributions, such as those from 
resonance decays and jet correlations. Those non-flow effects should not be included in the background to subtract 
from jet-like correlations. We shall refer to this jet-correlation background as flow-background. The anisotropic flow to 
be used for flow-background should be ideally that from two-particle method, vi{2~\ (ajf|, because jet-like correlation 
is analyzed by two-particle correlation method. Moreover, two-particle anisotropic flow contains fluctuations which 
should be included in jet-correlation flow-background 
. Non-flow is due to azimuthal correlations unrelated to the reaction plane, such as resonances, (mini)jets, or generally, 
£> ' clusters. Non-flow in two-particle t>2{2} was studied in a cluster approach [7] ; analytical form was derived for each non- 
flow component. In this paper, we shall demonstrate that the flow to be used in jet-correlation background subtraction 
1 should be the two-particle f2{2} excluding cluster correlations unrelated to the reaction plane, but including cross- 
terms between cluster correlations and cluster flow. We verify our result with Monte Carlo simulations. We then 
examine available jet-like correlation data from STAR [l[ and PHENIX Q with our refined flow-background. Finally 
we discuss implications of our result on the observed conical emission signal Q and on medium modification of jets 
in general. 



II. ELLIPTIC FLOW FOR JET-CORRELATION BACKGROUND 



In this section, we derive an analytical form for flow-background to jet-correlation in the cluster approach, as used 
in our non-flow study [Tj . We suppose a relativistic heavy-ion collision event is composed of hydrodynamic medium 
particles, jet-correlated particles, and particles correlated via clusters. Hydro-particles, high pt trigger particles, and 
clusters are distributed relative to reaction plane (ip) by 

^ = |^ [1 + 2« 2 cos 2(<£ -</>)] (1) 

with the corresponding elliptic flow parameter vi and multiplicity N. Particle azimuthal distribution with respect to 
a trigger particle is 



1 dN = dN hy dN d . k dN dtk dNddet 

N t dA<j) dA6 ^ dA<j> ^ dAcf) dA6 

k^jeteclus T kejet 



where Acf> = <j> — <f> t . In Eq. ([2]), dNdj e t/dA(j> is jet-correlation signal of interest. All other terms are backgrounds; 
Nd,k is the number of daughter particles in cluster k. 

If trigger particle multiplicity is Poisson and effects due to interplay between collision centrality selection (usually 
via multiplicity) and trigger bias are negligible, then the background event of the triggered (di-)jet should be identical 
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to any inclusive event, without requiring a high px trigger particle, but with all other event selection requirements as 
for the triggered event Q . Therefore, experimentally one can use inclusive events to obtain flow-background: 

1 dN _ ( dN hy i ^ dN^k v dNdA dN d , je t ( ] 



N t dA6 \ dAcj) ^ dAcj) ^ dA6 I dAcp 

k^jetddus k(£jet 

where a is a normalization factor, often determined by the assumption of ZYAM or ZYA1 (zero jet-correlated yield 
at minimum or at A<j) = 1) [Tol. [Tlj. and is approximately unity. 

In this paper we arc interested in the anisotropic flow to be used in jet-correlation background. So we will not 
concern ourselves with the background normalization, but only the background shape. We rewrite the background in 
Eq. © as 

dN bg _ dN hy dN d , k dN d , k _ dN hy y, dN d 

dA(j) ~ dA(f> ^ dAcj) ^ dAcj) ~ dAcj) ^ cl dAcj) ( ' 

k^jetZLclus k(£jet cl 

where we have eliminated subscript 'inc' to lighten notation. We have summed over all cluster types l cV including jet- 
correlation, where N c i is the number of clusters of type l cV . Different cluster types include jet and minijet correlations, 
resonance decays, etc. 

The hydro-background is simply 



dN hy = Nhy 

dAcj) 2tt 



(1 + 2v 2 ,tv 2 ,hy cos 2A0) (5) 



where v 2 .t is elliptic flow parameter of trigger particles and V 2 ,hy is that of hydro-medium particles. 
The cluster particles background is given by 

— — = / dcj> t pt((j>t) dcp k p cl (cf> k ) dA(p i f d {Acp i , c} k ) x — 6(A(pi + cj> k - Acj> - <j> t ) (6) 
dAcj) J J J 2tt 

where <f> t = <j> t - ip, 4>k = 4>k - ip, Afa = cj>i ~ cp k , and pt(<k) = ^ U + 2u 2 , t cos2^ t J and pd(4>k) = 

^1 + 2«2.c/ cos 2<J>kj are density profiles (i.e., i> 2 -modulated distributions) of trigger particles and clusters rela- 
tive to the reaction plane, respectively. We have assumed that the cluster axis (or cluster parent) distribution is also 
anisotropic with respect to the reaction plane. In Eq. ©, fd(A<j>i, cj> k ) = d J^ ( j ) k = is distribution of daughter 

particles in cluster relative to cluster axis (cluster correlation function) , which may depend on the cluster axis relative 
to the reaction plane cj> k [l2j]. Decomposing pt(4>t), we obtain 

dN 1 /'^TT plTT 2l) /•27T /"27T 

-jiri^Trl d(t>kPci(4>k) dA<f>ifd(dA<f>i,4> k ) + — — / dcj> k p cl (cj> k ) dAcj) i f d (Acj) u cj) k )cos2(Acj) i + cj> k -Acj)). 
dAcj) 2?r Jo J 2tt J J q 

( ? ) 

Because of symmetry fd(Acfii, cj> k ) = f d {-A(f>i,-<f> k ) and pd(4>k) = Pd(-<fik), we have 
f Q v dcf) k p c i(<t> k ) f Q v dAcj) l f d (Acj) l ,cj> k )sm2(Acj) l + cj> k ) = 0. Therefore 

— I dcj> k p c i(cf> k )Nd((j) k ) + — ^ cos2A0 / dcf>kPd(<f>k) dAfcf^Afc, (f> k ) cos2(A^ + <j> k ). (8) 



dAcj) 2tt J ™"-»vr*/ 2tt 
Realizing that elliptic flow parameter of particles from clusters is given by 

v 2 .d = (cos2(0 - ip)) cl = — / dcj> k p c i(cj> k ) dAcj) i N d (cj) k )fd(Acj) l ,4> k )cos2(Acj) i + cj> k ). (9) 

iV <i JO J0 

we rewrite Eq. ([5]) into 

</A,/ ^(l + 2 W2 , tW2 , d cos2A0). (10) 



dA<^ 2?r 

From Eq. ([5]) and (|10l) we obtain the total background as given by 



dl% = Nbg 
dAcp 2tt 



1 + 2w 2 ,i -r^-V2,hy + > 77 W 2jd COS 2A0 

V N b g <-f N h g I 



(11) 
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where 

N bg = N hy + J2 NdN d . (12) 

cl 

The v 2 s in Eqs. J3]), (fTU|) . and pT|) include fluctuations, so they should be replaced by yj {v\). The hydro-particles 
yj (v 2 ) is equivalent to two-particle «2{2} because there is no non-flow effect between hydro-particle pairs; same for 
the cluster yj (v 2 ) because there is no non-flow effect between different clusters (we consider sub-clusters to be part 
of their parent cluster). Thus Eq. (JTTJ) should be 

where 

Nhy rm , N c iN d 
V2.bg = -r F JL v 2 {2} hy + > — v 2 {2} d . (14) 

iVfeg ^ JN bg 

We note that here cluster includes single-particle (within a give px range) cluster, which generally is part of a parent 
cluster including particles of all px- Those single-particle clusters do not contribute to non-flow in v 2 {2}d, but they 
differ from single hydro-particles because they may possess different v 2 values. 

In principle, v 2j t should have a similar expression as Eq. (|14p out of symmetry reason: 

Nt,hy ro1 . N c l_tNt cl_t ro1 /-.r-x 

v 2 . t = 1 r^ JL v 2 {2} tMy + \ — - — ! v 2 {2} t , cl _ t . (15) 

cl_t 

where Nt,hy is number of high px trigger particles from hydro-medium (i.e., background trigger particles), t>2{2}t,hy 
is the elliptic anisotropy of those background trigger particles, N c i_ t is number of clusters of type ! cLi' containing 
at least one trigger particle, iV tjC /_ t is number of trigger particles per cluster, «2{2}t.ci_t is elliptic flow parameter 
of trigger particles from clusters, and N t ,tot = Nt,hy + / /Vd-t-Nt,d_t- The only difference is that trigger particles 

cljt 

are dominated by clusters (mostly jets), and those clusters are dominated by single-trigger-particle clusters; hydro- 
medium contribution to trigger particle population should be small. We note that jet-correlation functions are usually 
normalized by total number of trigger particles including those from hydro-medium background. 
If particle correlation in clusters does not vary with cluster axis relative to the reaction plane, 

v 2 {2} d = v 2 {2} c i(cos 2A0) cl , (16) 

and 

v 2 , bg = ^ 2 {2} hy + V ^^v 2 {2} cl (cos 2A0) d . (17) 

Nbg Nbg 



III. TWO-PARTICLE v 2 IN CLUSTER MODEL 



Obviously, the elliptic flow in Eq. (fl"4")) or (fl~7|) contains not only the two-particle anisotropy relative to the reaction 
plane, but also non-flow related to angular spread of clusters. How to obtain the elliptic flow as in Eq. (fT4|) or (|17|l ? 
In Q we have derived two-particle v 2 {2\ in a general hydro+cluster approach: 

vl{2} = (^v 2 {2U + Y.^ ± ^ 2 A +£^(<cos2A^) ci - U2 2 {2} d ). (18) 

The quantity in the first pair of parentheses in r.h.s. of Eq. (|18|) is elliptic flow due to correlation with respect to 
the reaction plane. The second term in the r.h.s. arises from cluster correlation Q; the small correction v|{2}(; is due 
to assumptions of Poisson statistics for number of clusters and particle multiplicity in clusters, but not the product 
of the two [7[. Since elliptic flow is formally defined to be relative to the reaction plane, the first term in r.h.s. of 
Eq. (fT5j) may be considered as "true" elliptic flow (except flow fluctuation effect), «a flow We note, however, it is not 
necessarily as same as hydro-flow because of contamination from clusters due to coupling between cluster correlation 
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and cluster flow. The second term in r.h.s. of Eq. (IT8l) can be considered as non-flow, f2. non-flow! non-flow is due to 
correlations between particles from the same dijet or the same cluster. Eq. (fT5|) can be expressed into 

w 2 {2} = ^2,flow + "2,noii-flow (19) 

Comparing Eq. (fT51) with Eq. (fT4")l . we see that 

V2.bg = Inflow, (20) 

i.e., the quantity in the first pair of parentheses in r.h.s. of Eq. (|18|) is the V2 parameter in Eq. (fT4|) that is needed in 
constructing jet-correlation background. In other words, elliptic flow parameter that should be used in jet-correlation 
flow background is the "true" two-particle elliptic flow (i.e., due to the reaction plane and including fluctuation). 



IV. MONTE CARLO CHECKS 



In this section, we verify our analytical result by Monte Carlo simulations. We generate events consisting of 
three components. One component is hydro-medium particles according to Eq. ([1]), given hydro-particles elliptic 
flow parameter V2,hy and Poisson distributed number of hydro-particles with average multiplicity Nhy The second 
component is clusters, given cluster elliptic flow parameter «2,ci and Poisson distributed number of clusters with 
average A^; each cluster is made of particles with Poisson multiplicity distribution with average Nd and Gaussian 
azimuth spread around cluster axis with ad- The third component is trigger particles with accompanying associated 
particles; the trigger particle multiplicity is Poisson with average Nt, and the elliptic flow parameter is V2,t- The 
associated particles are generated for each trigger particle by correlation function: 



/(A<£, <t>t) = Cfa)+ 5" S(< ^^ P 



N as (4> t ) 



2na as 



( 

exp 



V 



(A<f> - n + 6(4> t ) 



exp 



A^-7T-0(&) 



2al 



t) 



(21) 



where the near- and away-side associated particle multiplicities are Poisson with averages N ns (<pt) and N aa (^>t), 
respectively. The Gaussian widths of the near- and away-side peaks are fixed, and the two away-side symmetric peaks 
are set equal and their separation is fixed. All parameters in the jet-correlation function of Eq. (1211) can be dependent 
on the trigger particle azimuth relative to the reaction plane, <pt- 

We first verify Eq. (TT5)) by generating events with hydro-particles and jet-correlated particles. (We do not include 
other clusters except jet-correlations.) We use Nh y = 150, v 2 ,h y = 0.05, and N t — 2, w 2 ,t = 0.5. We use the large 
trigger particle V2 in order to maximize the effect of non-flow. For jet-correlation function, we generate back-to-back 
dijet with N ns = 0.7, N as = 1.2, a ns = 0.4, a as = 0.7, and 8 = (referred to as dijet model). We fix V2 in the 
simulation, i.e., V2 fluctuation is not included. We simulate 10 6 events and calculate w 2 {2} = (cos 2A^). Including 
only hydro-particles, we obtain t>2{2} ny = 0.05005 ± 0.00009, consistent with the input. Including all simulated 
events and all particles (hydro-particles and jet-correlated particles), we obtain U2{2}i nc = 0.05560 ± 0.00008. Using 
triggered events (events containing at least one trigger particle) only, we obtain t>2{2}trig_cvt = 0.05642 ± 0.00008. 
Using triggered events but excluding one dijet at a time (i.e., using the underlying background event of the dijet) 
and repeating over all dijets in the event, we obtain V2{2} Dg = 0.05568 ± 0.00008. We see that the background v 2 is 
as same as that obtained from inclusive events, «2{2}bg = V2{2}; nc , and both are smaller than that from triggered 
events only. 

We can in fact predict the inclusive event u 2 by Eq. (|17l) using the "hydro + dijet" model. The average yj (cos 2A0y) 



of jet-correlated particle pairs within the same dijet is \J (cos 2Axfiij) - et 
consistent with the expected value 



(cos 2A<f>) jet = 0.5054 ± 0.0004. This is 



(COS 2 A0) jet 



N ns + N a , 



■ exp 1 



-2al 



Nn 



N„ 



N a 



exp(-2cr 2 s ) cos 26 = 0.5046 



where 9 = 0. The average v /(cos2A0. 
equals to w 2 {2}d,jct = v2,t(cos2A0)jet 



■ij/ of pairs of particles from different dijets is 0.2516 ± 0.0008; 
0.5 x 0.5046 



it 



0.2523. The average (cos 2A(f)ij) for cross-talk 



pairs of background particle and jet-correlated particle is 0.10064 ± 0.00004; and it equals to the expected value 



5 



v/ «2{2}hy V2,t (cos 2 Ac/)) j C t = \/0.05 x 0.5 x 0.5046 = 0.1123. The inclusive event two-particle elliptic flow parameter 

is v 2 {2} = y/ (^jjjL x 0.05 + ^ x 0.2523) 2 + (0.5054 2 - 0.2523 2 ) = 0.05553; this is indeed consistent with 

v 2 {i\ mc or i>2{2}bg obtained from simulation. 

We now verify Eq. (|14l) or (ITT)) as the correct «2 to be used for jet-correlation background subtraction. We generate 
Poisson distributed hydro-particles with average multiplicity Nhy = 150 and fixed elliptic flow parameter V2,hy — 0.05. 
We generate Poisson distributed trigger particles with average trigger particle multiplicity N t — 2.0; we use different 
jet-correlation functions (discussed below). We also include clusters that do not have trigger particles (referred to 
as minijet clusters); the particle multiplicity per minijet cluster is Poisson distributed with average Nd = 5, and the 
number of minijet clusters is also Poisson distributed but we vary the average number of clusters Nd; we fix the 
cluster shape to be Gaussian with width ad = 0.5 (and average angular spread (cos2A0) c ; = exp (—2a- 2 .) = 0.6065), 
and also fix the cluster elliptic flow parameter v 2iC i = 0.20. We simulate 10 6 events and form raw correlation functions 
normalized by the number of trigger particles. In order to extract the real background v 2 from the simulations, 
we subtract the input jet-correlation function. If the jet-correlation function varies with the trigger particle angle 
relative to the reaction plane, the trigger multiplicity weighted average jet-correlation function is subtracted. We fit 
the resultant background function to B (1 + 2i>2,tV2,fit cos 2 A0) where B and v 2> at are fit parameters. We treat the 
input v-xt as known; we did not include any complication into i>2,t- We compare the fit V2 fit to the calculated one by 
Eq. ([TJ]) or (|17p . We study several cases with different shapes for jet-correlation function, as well as varying values 
for some of the input parameters: 

(i) "hydro + dijet" model: we generate back-to-back dijets accompanying trigger particles, without other clusters. 
The calculated v 2 , bg by Eq. is v 2 , bg = x 0.05 + x 0.5 x 0.5046 = 0.05500. 

(ii) "hydro + minijet + dijet" model: we include minijet clusters in addition to (i). The calculated v 2 ,bg by Eq. (fT7]) 
is v ^bg = acfl x °- 05 + Ml x °- 2 x °- 6065 + lirf x °- 5 x °- 5046 = 0.07126. 

(iii) "hydro + minijet + near-side + away-side double-peak" model: we generate jet-correlated particles by correlation 
function with double-peak away-side to replicate the experimentally measured reaction-plane averaged dihadron 
correlation function |l~3|, [l4| ■ We used the same Gaussian parameters for the correlation peaks as in (i) but 
= 1, thus (cos 2A</>) jet = 0.1689. The calculated v 2M by Eq. is v 2 . bg = ^ x 0.05 + x 0.2 x 0.6065 + 

x 0.5 x 0.1689 = 0.06813. 

(iv) "hydro + minijet + near-side + reaction-plane dependent away-side double-peak" model: we include reaction- 
plane dependent jet-correlation function similar to preliminary experimental data [l2[ • We have to use Eq. (fT4|) 
to calculate v 2 , bg , which gives v 2 , bg = x 0.05 + x 0.2 x 0.6065 + x 0.1423 = 0.06910. Note that, 
in this simulation of reaction-plane dependent jet-correlation signal, the number of jet-correlated particles is 
not 1.9, but rather 1.74. Also note that, due to the reaction-plane dependency of the jet-correlation signal, the 
elliptic anisotropy of jet-correlated particles cannot be factorized into the product of the trigger particle elliptic 
flow and the average angular spread of the jet-correlation signal as in Eq. [TH1 but has to be calculated by Eq. |H1 

We list our comparison in Table |U The fit i?2,fit is supposed to be the real background v 2tbg . The fit errors are due 
to statistical fluctuations in the simulation. As can be seen, the calculated v 2 , bg reproduces the real background v 2tbg 
in every case. The v 2jbg values differ from the hydro-background v 2 due to contributions from cross-talks between 
cluster correlation and cluster flow. Also shown in TableUare the two-particle «2{2} from all pairs in inclusive events. 
The i>2{2} values differ from v 2 ^ bg due to non-flow contributions between particles from the same dijet or the same 
cluster. 

Figure [Ha) shows the raw correlation function for case (iii) and flow background using the calculated v 2 ^ bg by 
Eq . [Til and normalized by ZYA1. Figure [TJb) shows the ZYAl-background subtracted jet-correlation function, using 
the calculated v 2 ^ bg for flow background. The background-subtracted jet-correlation is compared to the input signal. 
As shown in Fig.QJb), the shapes of the input signal and extracted signal are the same, which is not surprising because 
the calculated v 2 _ bg is the correct value to use in flow background subtraction. The roughly constant offset is due to 
ZYAl-normalization. 



V. HOW TO "MEASURE" JET-CORRELATION BACKGROUND (REACTION PLANE) v 2 

Two-particle angular correlation is analyzed by STAR and is decomposed into two components [TEj : one is the 
azimuth quadrupole, i>2{2D}, that is due to correlations of particles to a common source, the reaction plane; the 
other is minijet correlation that is due to angular correlation between particles from the same minijet or the same 
cluster. Of course, any such decomposition is model-dependent; because the functional form for minijet (or cluster) 
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TABLE I: Monte Carlo verification of analytical results of elliptic flow parameter to be used in jet-correlation background. 
Hydro-particle multiplicity, trigger particle multiplicity, jet-correlated near- and away-side multiplicities, number of minijet 
clusters, and particle multiplicity per minijet cluster are all generated with Poisson distributions, with averages Nhy, Nt, N ns , 
N as , N c i, and Nd, respectively. The jet-correlation function is given by Eq. <|2ip . with near- and away-side Gaussian width 
fixed to be a ns = 0.4 and a as = 0.7, respectively. The minijet cluster Gaussian width is fixed to ad = 0.5. The elliptic flow 
parameters for hydro-particles, trigger particles, and clusters are v 2t h y , v 2 ,t, and «2,ei, respectively, and are fixed over all events 
without fluctuation. We use Nh v = 150, Nt = 2, Nd = 5, v 2 ,hy = 0.05, and «2, c ; = 0.20. 



Case 



Parameters 



v 2 {2} 



V2,Rt 



Calculated v 2 



bg 



(i) hydro + dijet 

(ii) hydro + minijet + dijet 



(iii) hydro + minijet + near-side + 
away-side double-peak 



N cl = 0, v 2 , t = 0.5, 



1.2, 0.05557(8) 0.05505(8) 0.05500 




a ns = 0.4, a as = 0.7, 
N c i = 10, v 2 ,t = 0.5, 

C = 0, N ns = 0.7, N as = 1.2, 0.08465(6) 0.07115(8) 0.07126 
a ns = 0.4, a as = 0.7, 9 = 
N cl = 10, v 2 ,t = 0.5, 

C = 0, N ns = 0.7, N as = 1.2, 0.08172(6) 0.06815(8) 0.06813 
1 



(iv) hydro + minijet + near-side + 



<7„ s = 0.4, a as = 0.7, 
N c i = 10, v 2 ,t = 0.1, 

reaction-plane dependent away-side double-peak C = 0, N„ a = 0.7, N aa = 1.2, 0.08279 0.06883(35) 0.06910 
+ clusters a ns = 0.4, a as = 0.7, 6 — 1 




A<)) A<|) A()) 



FIG. 1: (a) Simulated raw correlation from the "hydro + minijet + near-side + away-side double-peak" model (Case iii in 
Table[T]). ZYAl-normalized flow-background using the calculated v 2 ,bg is shown as the curve. Note the input V2 in the simulation 
is purposely made much larger than real data to exaggerate non-flow effect, (b) Background-subtracted jet-correlation (data 
points) compared to the input correlation signal (upper curve). The background uses the calculated v 2 ,b g and is normalized 
to signal by ZYA1. The input signal shifted down by a constant is shown in the lower curve, (c) As same as (b) except the 
subtracted background uses the decomposed ti2{2D}. 



correlation is unknown a priori, one has to make assumptions about its functional form. Modulo this caveat, if cluster 
correlation and flow correlation is properly decomposed, the azimuth quadrupole should correspond to the first term 
in r.h.s. of Eq. (fl8)) . 

v 2 {2D} = ^v 2 {2} hy + £ ^v 2 {2} d . (22) 
N bg ^ JN bg 

This is identical to Eq. (|14[) . That is, the elliptic flow parameter from a proper 2D quadrupole-minijet decomposition 
is exactly what is needed for jet-correlation background calculation. 

A 2D quadrupole-minijet decomposition, with an assumption of the minijet correlation structure, has been carried 
out experimentally by STAR as a function of centrality but including all pt [Hj]. One may restrict to narrow pt 
windows to obtain «2{2D} as a function of pr, however, statistics can quickly run out with increasing pt because the 
2D decomposition method requires particle pairs. 

We perform a decomposition of flow and cluster correlation using our simulation data where the shape of cluster 
correlation is known from the simulation input. We form two-particle correlation between all particles (untriggered 
correlation). We fit the two-particle correlation with the sum of cluster correlation and flow to extract 1>2{2D} from the 
simulation data. Figure [He) shows the ZYAl-background subtracted jet-correlation function, using the decomposed 
t>2{2D} for flow background. The background-subtracted jet-correlation is compared to the input signal. The shapes 
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of the input signal and extracted signal are the same, which demonstrates that the decomposed i>2{2D} is close to 
the real elliptic flow value. Again, the roughly constant offset is due to ZYAl-normalization. 

One natural question to ask is why not to decompose jet-correlation and jet-background directly from high-py 
triggered correlation function. One obvious reason is, again, that jet-correlation shape is unknown a priori, thus 
one cannot simply fit triggered correlation to a given functional form. This is the same caveat mentioned above in 
particle pair correlation without a special trigger particle, where the minijet shape function has to be assumed in the 
decomposition of minijet correlation and flow. The situation in jet-like correlation is direr because the main interest 
of jet-like correlation studies is the investigation of medium modification to jet-correlation structure. Furthermore, 
even when the functional form of jet-correlation signal is known, as is the case in our simulation, we found that the 
decomposed jet-correlation signal shape deviates significantly from the input one. This is because the jet-correlation 
signal is not orthogonal to flow background, but rather entangled, both with near- and away-side peaks, and hence 
one can get false minimum \ 2 m decomposing the two components with limited statistics. 

VI. DISCUSSIONS AND SUMMARY 

In experimental analysis, vi values from various methods have been used for jet-correlation background. STAR used 
the average of the event plane ^{EP} and the four-particle «2{4} and used the range between them (or between V2{2} 
and «2{4}) as systematic uncertainties The event plane «2{EP} and two-particle vi\2\ contain significant 

non-flow contributions, while the non-flow contributions are significantly reduced in the four particle i>2{4} [16j. On 
the other hand, effect of flow fluctuation is positive in ^{EP} and «2{2} but is negative in «2{4}. This ensures that 
the true V2 is smaller than «2{EP} and «2{2}, and is most likely larger than t>2{4}. It is worth to note that the 
flow parameter to be used in jet-correlation background subtraction should include flow fluctuation effect as in «2{2}, 
which makes the «2{4} parameter as the lower limit rather conservative. 

The recently measured «2{2D} magnitudes from STAR are larger than «2{4} in peripheral and medium central 
collisions, confirming the validity to use «2{4} as the lower systematic limit of w 2 . In central collisions, however, the 
extracted f2{2D} is smaller than t>2{4} although the difference is significantly smaller than the difference between 
V2{EP} and ^{4}. This would suggest, assuming that the decomposed 1>2{2D} reflects the real flow background (i.e., 
the minijet shapes used in the decomposition is close to reality), that the used i>2 values for background calculation 
in dihadron correlation analysis in STAR would be too large by about ler systematic uncertainty. In three-particle 
correlation analysis Q, the U2{2D} was included in the vi systematic uncertainty assignment. 

STAR has also measured elliptic flow at mid-rapidity in the main Time Projection Chamber (TPC) using event- 
plane constructed by particles at forward and backward rapidities in the forward TPCs, V2{FTPC}. The obtained 
«2{FTPC} is smaller than ^{EP} using particles from the main TPC only, however, it is still significantly larger than 
«2{4}. This suggests that some but not all non-flow effects are removed from «2{FTPC}. The remaining non-flow may 
be dominated by the long range A77 correlation (ridge) observed in non-peripheral heavy-ion collisions [ll|, Ell Il8j - 
PHENIX used «2{BBC} results from the event plane method where the event plane is determined by particles in the 
Beam-Beam Counter several units of pseudo-rapidity away from particles used in jet-correlation analysis 19]. The 
rapidity gap in the PHENIX measurement of U2{BBC} is larger than that in the STAR measurement of ^{FTPC}, 
so non-flow effect should be smaller in the PHENIX measurement. However, it is possible that the W2{BBC} values 
used by PHENIX for background calculation can be also too large if non-flow ridge correlation persists to very large 
pseudo-rapidity gap. 

In summary, we have derived an analytical form for jet-correlation flow-background in a cluster approach. We 
argue that the elliptic flow vi parameter to be used in jet-correlation background is that from two-particle method 
excluding non-flow correlation unrelated to the reaction plane, but including cross-terms between cluster correlation 
and cluster flow. We have verified our result by Monte Carlo simulation for various jet-correlation signal shapes as well 
as varying other input parameters to the simulation. We demonstrate that the vi parameter to use in jet-correlation 
flow background is as same as the 1>2{2D} from a proper 2D quadrupole-minijet decomposition of two-particle angular 
correlation. However, we note that 2D quadrupole-minijet decomposition requires a model for minijet correlation 
shape, which gives rise to systematic uncertainty on the extracted t>2{2D} which require further studies. 
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